Proximity effect at the superconductor - topological insulator interface 
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We study the excitation spectrum of a topological insulator in contact with an .s-wave superconductor, starting 
from a microscopic model, and develop an effective low-energy model for the proximity effect. In the vicinity 
of the Dirac cone vertex, the effective model describing the states localized at the interface is well approximated 
by a model of Dirac electrons experiencing superconducting s-wave pairing. Away from the cone vertex, the 
induced pairing potential develops a p-wave component with a magnitude sensitive to the structure of the inter- 
face. Observing the induced .s-wave superconductivity may require tuning the chemical potential close to the 
Dirac point. Furthermore, we find that the proximity of the superconductor leads to a significant renormalization 
of the original parameters of the effective model describing the surface states of a topological insulator. 
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Topological insulators have emerged as an important class 
of materials characterized by robust quantum properties that 
can provide a significant potential for new device architec- 
tures. The recent observations of a two-dimensional (2D) 
topological insulator (TI) state in HgTe quantum wells HJ and 
of three-dimensional TI behavior in a family of bismuth-based 
materials, including Bi2Se3 and Bi2Te3 [2-6], represent a sig- 
nificant impetus to the experimental efforts in this area and 
a stimulus for the theoretical study of model systems with 
topologically ordered ground states ll7l- [T3l . The non-trivial 
topological properties of an insulator are intrinsically linked to 
the existence near the boundaries of gapless surface (or edge) 
states that are robust against disorder and interaction and have 
a characteristic Dirac type spectrum. 

Of special interest are the states that develop at the inter- 
face between a TI and an s-wave superconductor (SC), where 
the proximity effect generates a quasi-two-dimensional quan- 
tum state similar to a spinless p x + ip y superconductor, but 
without breaking time-reversal symmetry. This state supports 
Majorana bound states at vortices and was recently proposed 
by Fu and Kane [ 14] as a platform for fault-tolerant topologi- 
cal quantum computation. The proposal has generated a lot of 
interest due to the rather simple, weakly correlated nature of 
the underlying system. However, the existing calculations of 
the properties of this system are based on a phenomenologi- 
cal treatment of the proximity effect, where an s-wave pairing 
term is added to an effective Hamiltonian describing the TI 
surface. In addition, the general nature of the proximity effect 
between Dirac-type materials, such as graphene, and s-wave 
SCs has been a subject of great interest recently iflSHTTll . 

In this Rapid Communication we study the proximity effect 
at the interface between a TI and an s-wave SC starting from 
a microscopic model. In particular, we address and clarify the 
following questions: (1) how is the low-energy spectrum of 
the TI modified by the proximity effect induced at the inter- 
face? What is the dependence of this spectrum on the micro- 
scopic parameters of the system and what is the role of dy- 
namical effects? (2) What is the most general effective model 



for the interface and what is the dependence of its parameters 
on the relevant microscopic quantities? The answers to these 
questions have direct consequences for the stability of the Ma- 
jorana bound states that can be obtained at the interface. 

To study the proximity effect between an s-wave SC and 
the surface states of a strong TI, we consider a minimal mi- 
croscopic model defined by the following Hamiltonian: 

H to t - H-yi + Hsc + H t . (1) 

Here Hj\ is the TI term given by the tight-binding model on a 
diamond lattice with spin-orbit interactions 1 13 1, 

Hn = t J) c] a c jcr + U s o 2 • (djj X d;,)r:.r, ; . . (2) 

The first term is the nearest-neighbor hopping connecting the 
two fee sublattices of the diamond lattice. The second term 
connects second-order neighbors with a spin- and direction- 
dependent amplitudes. The direction dependence is given by 

the bond vectors d 1 ' 2 traversed between sites i and /'. The su- 

i j j 

perconducting term in Eq. ([TJ is given by mean-field Hamil- 
tonian 

H sc = ^ t^a^a^ + ^(Aoa^ + h.c.) (3) 

(if), a i 

The model is defined on a simple hexagonal lattice with a lat- 
tice constant for the triangular basis that ensures simple in- 
terface matching conditions between the SC and the (1,1,1) 
surface of the TI. In general, the nearest neighbor hopping t' tj 
has different values along directions parallel and perpendicu- 
lar to the interface. The tunneling term in Eq. ([T]) is 

ff t = 2 ^ a L c .^ + 4^), ( 4 ) 

(0) 

were t is the tunneling matrix element that characterizes the 
transparency of the interface between TI and SC. The fermion 
operators cj lT and a llT operate in the Hilbert space of the TI and 
superconductor, respectively. 
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FIG. 1: (Color online) Band structure for a slab with a (1, 1, 1) sur- 
face of a model TI described by Eq. {2|. The blue (dark gray) areas 
represent bulk states and the in-gap lines are surface states. The right 
panel shows a zoom-in of the anisotropic Dirac cone at the M point 
of the surface Brillouin zone. 

In the absence of tunneling, the spectrum of the model 
TI described by Eq. Q develops a bulk gap, provided an 
anisotropy is introduced in the hopping along the four near- 
est neighbor bonds, t — > t + 6t a , a — 1,4 0~3). In the presence 
of a boundary, gapless surface sates develop inside the bulk 
gap. The results for a slab with a (1, 1, 1) surface and param- 
eters t = 1, Aso = 0.3 and St = 0.25 for bonds with a surface 
projection oriented along M (which we choose as the x di- 
rection) are shown in Fig. [T] The surface states spectrum is 
characterized by presence of one (anisotropic) Dirac cone in 
the vicinity of the M point, which represents the signature of 
a strong topological insulator ifTSl . 

Next, we discuss the effective model that describes the TI 
surface states, as a preliminary step for identifying an effective 
model for the TI-SC interface. Notice that, for each wave vec- 
tor k in the vicinity of M there are exactly two surface states 
with energies si(fc) = e^(fc) and sjik) = e N+1 (k), where e n (k) 
with n = 1, . . . , 2N are the energy eigenvalues for a slab con- 
sisting of N layers. The corresponding wave functions are 
U\,k( r , c) an d u 2,k( r ' °~)> respectively. An effective Hamilto- 
nian describing the surface states can be obtained by project- 
ing the full TI Hamiltonian (|2]) onto the subspace spanned by 
u\ t k and «2,fc. If we choose a basis for this subspace defined 
by certain wave functions u ± k, an arbitrary surface state can 
be expressed as 

<Afc(T, cr) = ^ fA,kUA,k( r > °")> ( 5 ) 

,i=± 

where ip ± ^ depend on the two-dimensional momentum k. Af- 
ter projection, the Schrodinger equation for the surface states 
becomes an equation for the spinor (p(k) = (ip+j~,<P-,k) T or, 
in real space, for its Fourier transform <p(r), where r is now a 
2D vector. Thus, the full three-dimensional quantum problem 
for the TI surface states is reduced to an effective 2D problem 
defined by H eS (k; A, A') = {u, lk \H T i\u A ^ k ). 

However, choosing an appropriate basis for this subspace 
has to take into account two requirements that stem from con- 
siderations of symmetry and analyticity. First, it is convenient 
to have an effective Hamiltonian with manifest time-reversal 
(TR) properties. To this end, one can choose a pseudo-spin 



basis « ±j fc(r, cr) with the property &u ± k = ±M+,fc, where 
= ia-yK is TR transformation operator and K denotes 
the complex conjugation. Under such a choice, the spinor 
<p(k) = (<p + k,<f-,k) T satisfies the relation ®<p(k) = icr y (f)*(-k), 
which is identical to the time-reversal transformation of phys- 
ical spinors. Given that spin is not a good quantum number in 
the presence of spin-orbit interactions, it is natural to construct 
the effective Hamiltonian using the pseudo-spin basis u^^m- 
Note that the eigenfunctions U\k and u^k have complicated 
TR properties, hence they do not represent a convenient basis. 
Second, we may be interested in a real space formulation of 
the effective theory that would apply to spatially inhomoge- 
neous situations, for example, when impurities or vortices are 
present. Thus, it is desirable to construct an effective model 
using real space spinors (p(r) = {ip + {r), tp_{r)) T . A surface 
state can be expressed in terms of Wannier functions as 

«A(r, cr) = J] ^(R)«u(r - R, cr), (6) 

,t=±;R 

where w^r - R, cr) = £ k e' k (r_R) « /1 k (r, cr) and u ±y k are appro- 
priate combinations of the energy eigenstates ui^-k, so that 
they transform as pseudo-spins under time reversal. To con- 
struct exponentially localized Wannier functions |T8l . the ba- 
sis states u ± k(r, cr) have to be analytic in the vicinity of the M 
point where energy levels cross. Note that the limit k — » M is 
ill-defined for the energy eigenstates u, h k, as it depends on the 
direction along which the M point is approached. 

A simple way of constructing the pseudo-spin basis actually 
used in our numerical analysis is to form linear combinations 
of ui,2;k( r , 0") that are spin-polarized in the top TI layer, 

(u + , k (r N ,1)\ ff .M\ (u-,k(rN,V\ t( .x/0\ _ 

U,*(r„,l)] = /( >)' 1^4)J =/( - fe VJ' () 

where /(fc) is a positive smooth function. Generally, the func- 
tions u±,k(fcr) are not orthogonal, but one can easily construct 
almost spin-polarized orthogonal wave functions as 

«+,fc = N+(fc) (S +>fe + a(k)u- h ) , (8) 
u-k = N-(k) (S_ k + a*(k)u+,k) (9) 

where a(k) = V 1 - K«+,fel«-,fc)l 2 ) /<«+,fel"-fe> van- 

ishes at the Dirac point and N ± (k) are normalization factors 
chosen so that u + j,(r N ,1) and U-j^TnA) be positive. The 
functions u ± k obtained using this prescription have the trans- 
formation properties of spinors under time reversal and are 
smooth in k around the M point. The construction is valid for 
all values of k corresponding to surface modes. 

If we are only interested in the properties of the surface 
states in the vicinity of the M point, the effective 2D model 
can be constructed using the pseudo-spin basis at the M point, 
u±,k M > an d k ' P perturbation theory. One obtains 

Hfj ^E M + Bik) ■ cr, (10) 

where E M is the energy at the M point, B(k) is a linear func- 
tion of k and cr = (cr x , o- y , cr z ) are Pauli matrices represent- 
ing the pseudo-spin. The most general form for the effective 
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FIG. 2: (Color online) Spectrum of the full Hamiltonian {TJ| for a 
TI slab with (1, 1, 1) surfaces having an interface with an s-wave 
superconductor. Bulk TI states are green (light gray) and SC states 
are blue (dark gray areas). The proximity effect induces a gap in the 
surface states spectrum at the interface (red/dark gray in-gap lines), 
while the states on the opposite surface have an unperturbed Dirac 
cone dispersion (orange/light gray lines). The light blue (light gray) 
in-gap points are calculated using an effective theory for the interface 
described by Eqs. (JTTJ, and l |14|16[ l. 

Hamiltonian can be determined using the time-reversal prop- 
erties of the pseudo-spin basis u ± ^,'- 



Hfj = L (k) + J] L,(k)cr ; - 



(11) 



i=l,2,3 



where Lo(k) and L,(k) are even and odd smooth functions of 
k, respectively. Unlike, the Hamiltonian ( [T0| derived from 
k • p perturbation theory, the above Hamiltonian is valid for 
all momenta over which the surface states are defined. For a 
generic anisotropic TI model an important feature is the pres- 
ence L z (k), which leads to unequal spin population. Note that 
in the limit of k — > M we recover the Hamiltonian ( fT0| ). 

Next, we turn on the coupling between the TI and the SC. 
As a result of the proximity effect, a finite gap opens in the 
surface states spectrum. The result of the numerical calcu- 
lation for the full tight-binding Hamiltonian ([T]l is shown in 
Fig. [2] for the same TI model parameters as in Fig. [T[ a 
superconductor with fjj = 0.5f and t' ± = 0.45 for hoppings 
parallel and perpendicular to the interface, respectively, and 
Ao = O.lf and an interface characterized by a transparency 
t = 0.3f. To better understand the dependence of the proxim- 



ity effect on various model parameters, it is useful to identify 
the SC contribution to the interface effective model and add it 
to TI contribution given by Eq. ( fTT| ). Within the mean-field 
model given by Eq. the superconducting degrees of free- 
dom can be integrated out and replaced by the interface self- 
energy Y. aa -(r,r' ,u>) = -T(r,ri)G CT a-irur 2 ;aj)T t (r2,r'), 
where Ga-a-'ifi, r^, co) is the superconducting Green's func- 
tion in the Nambu space and f is the matrix describing the 
tunneling between top layer of the TI (layer N of the TI-SC 
heterostructure) and bottom layer of the superconductor (layer 
N+l). For an s-wave superconductor, we have 



E(fc,w) 



LOT {) + 6„ ,(fc)T z + A T X 

e m (ky 



+ A 2 Q - a) 1 



\Xmk(ZN+l)\ (12) 



with Xmk(zN+\) being the value on the bottom layer of the 
energy eigenstates in the metallic phase. This self-energy in 
terms of the partial density of states at the interface, v(s, k) = 
5{E-e m (k))\x m k{zN+\)^ \ which in our model takes the sim- 
ple form v(s, k) = ^ Jl - l£= fr^, where A = 2t x is the 
bandwidth in the metallic phase of a slab described by the 
Hamiltonian ([3]) and e(fe) is the energy value at the middle of 
the band for a given momentum parallel to the interface. For 
surface states in the vicinity of k m k^ and assuming Ao <s A, 
we have 



S(fe,£U) 



(U>T + A T A ) 



(13) 



where A = 2f/A^l - (e(k M ) - fj) 2 /A 2 , with p. being the 
chemical potential of the metal. Note that this approximation 
ignores the momentum dependence of the self-energy induced 
by the s-wave superconductor. 

To construct the effective model for the proximity effect, we 
project the above equation onto pseudo-spin basis M±,fc(r,cr). 
The effective self-energy at the interface becomes 

FfJ(k, oS) = ^ u Ak (r N , cr)S^,(fe, (i))u x -jt(r N , cr') 



Aco 



XI U *k( r *f' cr ) M -l',fc( r A'' °") (I 4 ) 



F ( ^,(k,to) = X Mafe(rj V ,cr)E^,(fe,w)M / i',-fe(r w ,o-') 



AAn 



V^ 7 



^ UAk(r N , cr)uv-h(r N , -cr). (15) 



Here T, N ^ A are the normal (diagonal) and anomalous (off- 
diagonal) parts of the self-energy and k = (k x , k y ). In general, 
both Fx9, and F\), have even and odd components in k. The 
anomalous part of F\ x (k, to) is an even function of frequency 
and has singlet A,(k, to) = A s (— k, to) and triplet A,(k, to) - 
-A t (— k, cS) components. Note that the relevant fc-dependent 
corrections arise from the pseudo-spin eigenstates u A j-(r N , cr). 
Even though the construction of the pseudo-spin basis is de- 
fined up to a ^-dependent unitary transformation that preserve 
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FIG. 3: (Color online) Surface state localized at the interface between 
a TI and an s-wave SC. For zero tunneling (not shown), the state 
resides entirely inside the TI and decays exponentially away from 
the boundary (up to an even-odd oscillatory factor). At finite trans- 
parency, the surface state penetrates inside the SC. The spin-down 
(lower amplitude) component vanishes on the top TI layer, thus it 
does not propagate inside the SC if only nearest-neighbor tunneling 
is considered. 

the time-reversal properties of the pseudo-spin basis, such 
transformation preserves |A s (k, to)\ 2 and |A,(k, to)\ 2 . We em- 
phasize that the emergence of the triplet component can also 
be viewed as a property of the k • p perturbation theory for 
the Bogoliubov-de Gennes (BdG) Hamiltonian. This property 
is generic and independent of the position of the Dirac point. 
Finally, the spectrum of excitations at the interface can be de- 
termined by solving the BdG equation 



Det(i/^'(fc) + F(k, a>) - u>) = 0. 



(16) 



A comparison between the spectrum of the TI-SC het- 
erostructure described by Eq. ([TJ and the spectrum of the ef- 
fective model in Eq. ( 16 1 is shown in Fig. [2]for a system with 
a chemical potential that crosses the Dirac point, fi — 0. Sim- 
ilar calculations were performed for systems with /j + and 
for systems with a different interface orientation, i.e., (100) in- 
stead of (1 1 1). Adding a potential barrier at the interface was 
also considered but the results remain qualitatively the same. 

We now discuss the main implications of the microscopic 
calculations described above. In the limit k as Icm and y m 
the triplet component A, vanishes due to the time-reversal 
symmetry and pairing becomes purely s-wave, as conjectured 
in Ref. Ifl4ll . The low-energy spectrum can be obtained by 
expanding the self-energy to lowest order in k and to 



Det 



1 



l+A/A 



(ucrk)T z + A 



A + A 



0. (17) 



We find that in this limit, consistent with previous phe- 
nomenological results Ifl9l , the form of the above effective 
BdG equation resembles the equation for the proximity ef- 
fect at the interface of a TI and an s-wave superconductor |[T4l 



with the critical difference that the effective parameters of the 
TI surface states such as the velocity v and chemical potential 
fi are strongly renormalized by l+A/A. This renormalization, 
which can be thought of as a consequence of the penetration 
of the TI wave function into the SC (see Fig. H), determines 
the energy scale for all subgap excitations at the interface, e.g., 
the Majorana modes. In particular, it determines the size of the 
minigap that protects the Majorana bound state, hence it has 
direct consequences for the robustness of topological quantum 
computation 1 19]. 

As the chemical potential is tuned away from the Dirac 
point fi m the triplet component of the pairing potential de- 
velops. This follows from the fact that the wave functions 
on the TI are not spin-polarized due to the spin-orbit cou- 
pling as shown in Fig. [3] The nature of the pairing away 
from the Dirac point can be analyzed by considering the BdG 
Eq. (16 1, for to, A «: Ao (to ignore retardation effects for 
simplicity), and keeping the lowest-order k dependence as 
Det([v,- • ko-j-fj]T z + /l[l+u,' • ko-j]T x ) = 0. Thus in addition 
to the proximity-induced singlet pairing of magnitude A, there 
is a proximity induced spin-triplet pairing of magnitude A\u\kf 
where kp~y.lv is the Fermi wave-vector on the Dirac cone. 
This triplet pairing can become significant when fi > uA/u 
and can lead to a closing of the superconducting gap on the 
surface state of the TI, which has important consequences for 
topological quantum computation. In particular, the resulting 
proximity induced superconducting state can be qualitatively 
different from the one where the chemical potential is near the 
Dirac point and is not expected to support Majorana fermions. 

We calculated the excitation spectrum of a topological insu- 
lator in contact with an s-wave superconductor, starting from a 
microscopic model, and we developed an effective low-energy 
model for the proximity effect. We found that a /?-wave pair- 
ing component generally leads to the reduction in the induced 
quasiparticle gap. We also showed that virtual propagation 
of the electrons into the superconductor leads to a signifi- 
cant renormalization of the effective surface TI Hamiltonian. 
These findings have important implications for topological 
quantum computation. 

This work is supported by DARPA-QuEST, JQI-PFC and 
LPS-NSA. 
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